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^ , Abstract 

We study the time evolution in system of N bosons with a relativistic dispersion law interacting 
through a Newtonian gravitational potential with coupling constant G. We consider the mean 
2 ' field scaling where N tends to infinity, G tends to zero and A — GN remains fixed. We investigate 

the relation between the many body quantum dynamics governed by the Schrodinger equation and 
the effective evolution described by a (semi-relativistic) Hartree equation. In particular, we are 
interested in the super-critical regime of large A (the sub-critical case has been studied in [^[T^). 
^ ' where the nonlinear Hartree equation is known to have solutions which blow up in finite time. To 

lO ! inspect this regime, we need to regularize the interaction in the many body Hamiltonian with an 

CO . N dependent cutoff that vanishes in the limit N — >■ oo. We show, first, that if the solution of the 

nonlinear equation does not blow up in the time interval [— T, T], then the many body Schrodinger 
dynamics (on the level of the reduced density matrices) can be approximated by the nonlinear 

|0 I Hartree dynamics, just as in the sub-critical regime. Moreover, we prove that if the solution of the 

nonlinear Hartree equation blows up at time T (in the sense that the H^^^ norm of the solution 
diverges as time approaches T) , then also the solution of the linear Schrodinger equation collapses 
(in the sense that the kinetic energy per particle diverges) if t ^- T and, simultaneously, N —>■ oo 
sufficiently fast. This gives the first dynamical description of the phenomenon of gravitational 

^ ■ collapse as observed directly on the many body level. 

^: 

1 Introduction and main results 

We consider systems of gravitating bosons known as boson stars. Assuming the particles to have a 
relativistic dispersion, but the interaction to be treated classically (Newtonian gravity), we arrive at 
the A^-particle Hamiltonian 
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acting on the Hilbert space L^(M^^), the subspace of L'^{M.^^) containing all functions symmetric 
with respect to arbitrary permutations (here we use units with h = 1, c = 1, and m = 1, where m 
denotes the mass of the bosons). 



We are interested in the mean field limit where A^ — )• oo, G — )• so that NG =: A remains fixed. 
In other words, we are going to study a family of systems, parametrized by the number of bosons 
N, described by the A^ particle Hamiltonian 



F^. = V./i-A.,-Av^^. fi.ii 



N , N 

^^ V ^ N '^-^ \xi — Xj\ 

The system is critical, and it behaves very differently depending on the value of the coupling 
constant A > 0. The criticality of the system is a consequence of the fact that the kinetic energy 
scales, for large momenta, like the potential energy (both scales as an inverse length). The potential 
energy can be made arbitrarily large (and negative) by moving the particles closer and closer together 
(iV particle in a box of volume i^ have a potential energy of the order N£~^, taking also into account 
the 1/A^ factor in front of the interaction energy). However, in order to localize particles in a small 
volume we have to pay a price in terms of kinetic energy (to localize A'^ particles within a box of 
volume i^, we need an energy proportional to Ni~^). This simple observation implies that, for 
small values of the coupling constant A, the kinetic energy dominates the potential energy, and 
that, for sufficiently large A, the kinetic energy needed to bring particles together is not sufficient to 
compensate for the gain in the potential energy. 

For every N £ N, there exists therefore a critical coupling constant Acrit(A^) such that H]\f is 
bounded below for all A < Acrit(A^) and such that 

mf — - — — — = — oo 

^eL2(M3iV) 



for all A > Acrit(A^)- It was proven in [T7] that the critical constant is given, as A^ — >• oo, by the 
critical coupling constant for the Hartree energy functional 

£u.nrM = / d. (1 - A) VV(.) '-U dxd.M^l^^ . (1.2) 

J 2 7 \x-y\ 

More precisely, it was proven in p/7j that, as A^ — )• oo, Acrit(A^) — ^ ^Sif where 

1 _ lfdxdy\ip{x)\'^\ip{y)\'^\x-y\-^ 



Aj^rit <peL2(R3),||^||=i 2 /dx||V|V2y,(j;)|2 

Note that, with this definition, SnuTtreei^) > for ah y? G H^/'^{R^) if A < A^^^ while 

inf f Hartree (V^) = "OO 

if A > A^-^. It is also possible (see [I^) to give bounds on the fluctuations of Acrit(A^) around A^^^.: 

A«it(l - ciAr-V3) < x^^.,,(N) < A«it(l + C2N-') 

for appropriate constants ci,C2 > 0. The value of A^^^. is not explicitly known. By Kato's inequality, 
\xi — Xj\^^ < (7r/2)|V^.|, it is easy to see that A^j|. > (4/7r) ~ 1.3. In |l8l[T7], it is also shown that 
A^rit < 2.7. 

Let us discuss first the subcritical case A < A^j^. In this case, the Hamiltonian p.ip has, at 
least for sufficiently large A^ (so that A < Acrit(A^))) a unique realization as a self-adjoint operator on 



L^(R^ ) and therefore generates the one-parameter group of unitary transformation Uj\f{t) = e ^^^t ^ 
with t G M. The unique global solution of the A^-particle Schrodinger equation 

idt^N,t = Hn^n,U i^N,t=0 = V'iV G L2(R3^) (1.3) 

which governs the time evolution of an arbitrary initial A^-particle wave function ip^- is then given 

by 1pN,t = UNit)lpN- 

Consider now the time evolution (jl.3|) of a factorized initial data tpN' = ^^^ for some (p G L^(M'^) 
(here we use the notation c/?®^(x) = Y[j=i 'Pi^j)^ where x = {xi, . . . , xn))- Of course, factorization is 
not preserved by the evolution. Nevertheless, because of the mean field character of the interaction, 
one may expect factorization of the evolved wave-function ipN,t to be restored, in an appropriate 
sense, in the limit of large A^. If we assume, formally, that 

i^N,t ^ ^r , (1.4) 

then it is easy to derive a self-consistent equation for the evolution of the one-particle orbital <pt- 
In fact, if (jl.4p is correct, the potential experienced by the particles can be approximated by an 
averaged, mean field, potential given by the convolution —A |.|~^ * Iv'tP- Therefore, ()1.4p implies that 
ipt must evolve according to the semirelativistic nonlinear Hartree equation 

idt<pt = \/l-A(^t-A( — * IvJtp ) ipt (1.5) 



with initial data (pt=o = f- The question now is in which sense can the factorization (11. 4p hold true. 
It turns out that (jl.4p should be understood on the level of the marginal densities associated with 

tpN,t- 

Let 7Ar^i = \ipN,t){ipN,t\ be the orthogonal projection onto the solution of the A^-particle Scrodinger 
equation ipN,t = e~*^^*^7v with factorized initial data tpN = ^^^ ■ Then, for A: = 1, . . . , A^, we define 
the A:-particle marginal (or reduced) density jj^^ associated with ^7v,t by taking the partial trace of 
lN,t over the last N — k particles; that is 

(k) 
lN,t = T^^k+l,k+2,...,N lN,t ■ 

In other words, 7)^^ is defined as a non-negative trace class operator on L'^(M. ) with kernel 

7^,i(xfe,yfc) = / d^N-k7N,t{^k,XN-k;yk,^N-k) 

dxAr_fc^Af,f(xfc,XAr_fc)V'Ar,t(yfc,XAr_fc), 



where we introduced the notation x^ = (xi, . . . , Xk), Yk = {iJi, ■ ■ ■ , Uk), ^N-k = {^k+i, ■ ■ ■ , xn)- 

The first rigorous proof of the validity of (jl.4p on the level of the marginal densities (and hence 
the first derivation of ()1.5p ) has been obtained, for the subcritical regime, in [2]. More precisely, 
under the condition that A < A/n (which is smaller than A^^^), it is proven in [2] that, for every fixed 
A; > 1, and for every t G M, 

Tr ^(J)_|^^^^^^|0fc ^0 (1.6) 



as A^ — )• oo. Here ipt is the solution of the semirelativistic nonhnear Hartree equation (jl.Sp with initial 
data </3t=o = ^- Note that the convergence of the marginal densities implies that, for an arbitrary 
/c-particle observable J^^\ we have 



(V'iv,, (J('') ® i(^-'=))V'^,> ^ (^f , j('Vf 



k\ 



as A^ — 7- oo. In this sense, the solution of the A^-particle Schrodinger equation ipN,t can be approxi- 
mated, for large N , by products of the solution of the one-particle semirelativistic Hartree equation 
(jl.Sp . Observe here that the semirelativistic Hartree equation (jl.Sp is locally well-posed in the en- 
ergy space /7^'^(R^), for arbitrary A S M. In fact, it is proven in |16] that, for every c^ G H^/'^iM?), 
there exists a maximal < T < oo and a unique solution 93. G C((— T, T),i^^'^(IR^)) of (II. Sh in the 
time interval t G {—T,T). Here, either T = 00 (and then the solution is global), or T < 00 and 
11^5411/^1/2 — )• 00 as t — )• r or as t — )• —T (in this case, ipt blows up in finite time). For A < A^j^, the 
unique solution ipt is also shown to be global (hence T = 00 in this case); for this reason, (jl.6p makes 
sense for all t G M. 



Note that the result (jl.6p is just one of the several results concerning the derivation of effective 
evolution equations from first principle quantum dynamics which have been obtained in the last 
years. Other results of this type concern the derivation of the non-relativistic Hartree equation 
in the mean-field hmit (see for example [22l El H EU HOl HSl dU IS]) and the derivation of the 
(non-relativistic) Gross-Pitaevskii equation for the description of the dynamics of initially trapped 
Bose-Einstein condensates (see [U [5l El [7] and, very recently, |20]). Although most of these papers 
deal with non-relativistic particles, the authors of |19] also consider semirelativistic bosons interacting 
through a Newtonian potential, in the subcritical regime A < A^j(.. They improve ()1.6p . by giving 
an explicit bound on the rate of the convergence; more precisely, they show that 



Tr 



7S-I^t)(9^*r <% (1.7) 
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for a constant Ck,t = Ct Vk, where Ct grows at most exponentially in t G M (under additional 
assumptions on the dispersion of ipt, Ct is bounded uniformly in t). 

So far, we considered the subcritical regime A < A^j^. Let us discuss now the supercritical regime 
A > A^j^. The criticality of the Hartree energy functional remarked earlier can also be observed on 
the level of the time-dependent semirelativistic nonlinear Hartree equation (|1.5p . On the one hand, 
(II. Sp is globally well-posed for A < A^;^. On the other hand, it turns out that, for A > AJj^^^., (II. 5p has 
solutions that blowup in finite time. More precisely, it was proven in [9j that, for every spherically 
symmetric Lp G H^{^) with <?Hartree('/') < such that |||2;|99|| < 00, the unique maximal solution 
ipt G C{{—T,T),H^''^{E?)) of (jl.Sp with initial data 994=0 = 'f blows up in finite time, in the sense 
that T < 00 and 

llv^'i 11/^1/2 —5- 00, as t — ;■ T orast— )•— T. 

Solutions of ()1.5p exhibiting blowup in finite time are supposed to describe, within the framework 
of Chandrasekhar's theory, the gravitational collapse of bosons stars. The expectation that blowup 
solutions of the semirelativistic Hartree equation describe the collapse of boson stars is based on the 
unproven assumption that the many-body dynamics can be approximated by the Hartree dynamics 
also in the supercritical regime A > A^j^. and all the way up to the time of the (nonlinear) blowup. 
In this paper, we give a rigorous proof of this physical assumption. We show, first of all, that the 
convergence ()1.7p also holds in the supercritical case A > A^j^, if the norm ||99s||^i/2 stays bounded in 
the interval [0, t]. Moreover, we prove that the convergence of the A'^-particle Schrodinger evolution 



towards the Hartree dynamics does not only hold in the sense of ()1.7p : instead, it also holds (again 
assuming that the norm ||c/3s||//i/2 remains bounded in [0, t]) with respect to the stronger energy norm, 
at least for the one-particle marginal density (this means that (1 — A)^' ^^^^{1 — A)^'^ converges 
to (1 — A)-^'^|(/9t)((/9t|(l — A)^'^ in the trace norm, as A^ — )• oo). Note that this is the first proof of 
the convergence of the many body Schrodinger evolution towards the Hartree dynamics with respect 
to the energy norm, not only for supercritical semirelativistic boson stars, but for any mean field 
system. As a consequence of the convergence in energy, we show that the solution of the Hartree 
equation (|1.5p really describes the collapse of the many body system. 



Let us now describe our results in more details. We are interested in the dynamics generated 
by the iV-particle Hamiltonian (jl.ip in the supercritical regime A > A^^^. (although at the end our 
results will also hold in the sub-critical regime). The first issue that we have to face is that, since the 
form defined by H^ is not bounded below, the Hamiltonian H]\f does not necessarily have a unique 
realization as a self-adjoint operator on the Hilbert space L^(M^^) (and thus the time-evolution is 
not necessarily well-defined). For this reason we choose a sequence a = {aN)N>i with a^ > for 
all A^ G N and ajsf — t- as A^ — )• oo, and we define the regularized A^-particle Hamiltonian 

H% = y,fi^A^.--y, — ^ — . (1.8) 

^ Z^V ^' N ^ Xi-Xi +aN 

The regularized Hamiltonian H'^ defines now a quadratic form on Lg(M'^ ) which is clearly bounded 
below, for every A^ G N, since 

By Friedrichs theorem, H"^ has a unique extension as a self-adjoint operator on L^(M^'^), with domain 
H^/'^i^^). Hence H'^ generates the one-parameter group of unitary transformations U^{t) = 
g-itH"^ ^ f ^]^^ and therefore the A^-particles Schrodinger equation 

idttpN,t = H'^ijJN,t with initial condition ipN^t=o = '4'N (1-9) 

is globally well-posed (it has the unique solution V'tv,* = e~'^^''^^ipN^ for all t G M). 

From the physical point of view, the introduction of the cutoff a is justified by the observation 
that on very short length scales, the Newtonian potential is effectively regularized by the presence of 
other forces (such as electromagnetic or nuclear forces) or because of general relativity effects. The 
results that we will state and prove below concern the limit of large A^ and small a^. How fast ajv 
tends to zero is irrelevant to establish the convergence of the Schrodinger evolution to the Hartree 
dynamics in the trace norm, analogously to (|1.6p (although, of course, the rate of the convergence 
depends on a^)- On the other hand, to show convergence in the energy norm, we will need to assume 
that aTv does not converge to zero too fast; more precisely, we will suppose that there exists /3 > 
such that N^a^ — )• oo. This condition, which allow for any power law decay, still leaves a lot of 
freedom in the choice of a^v (physically, conditions on the decay of a translate into restrictions of the 
range of systems for which the approximation of the many body evolution by the Hartree dynamics 
is applicable). 

We study the time evolution generated by the regularized Hamiltonian (II. Sp on factorized initial 
data t/jN = if® for ip G -ff^(M^). We compare the marginal densities associated with the solution of 
the A^ particle Schrodinger equation il)N,t = e~*^^*'(/'A' with products of the solution to the Hartree 
equation (II. 5p . Note that the cutoff disappears in the limiting Hartree equation, because of the 
assumption that a^ — )• as A^ — )• oo (part of the proof of the convergence will consists in estimating 



the distance between the solution ipt of (jl.Sp and the solution ipf' of a regularized Hartree equation 
with interaction — A/|j;| replaced by the regularized interaction — A/(|j;| + a) in the limit of small a). 

The first main result of this paper is the following theorem. Under the assumption that the 
solution ift of ()1.5p has a bounded i^^' ■^-norm in the interval [— T, T] (which means that there is no 
blowup, up to time T), we prove the convergence of the marginal densities associated with the solution 
ipjy^t of the A^-particle Schrodinger equation (jl.9p to the orthogonal projections onto products of (ft- 
The theorem also gives an explicit upper bound on the fluctuations around the Hartree dynamics. 

Theorem 1.1. Fix X € R, if e H'^{R^) with \\ip\\ = 1 and set tpN = 93®^- Let iljN,t = e-^^^^Af ^e 
the evolution of the the initial wave function tp^ with respect to the Hamiltonian /il.8]) . and let '^\q ^ 
he the one-particle reduced density associated with ipN,t- 

Denote by ipt the solution of the nonlinear Hartree equation M.5\) with initial data (pt=o = ^- F^^ 
T > such that 

K := sup ||i/3j||j:^i/2 < oo . (1-10) 

\t\<T 

Then there exists a constant C = C{k,T, W^Wh^) such that 

Tr|7«-|(pt)((^t|| < c[-j=+aM). (1.11) 




for all t €z M with \t\ < T, and for all N sufficiently large. In particular, if on — )• as A^ — )■ oo, it 
follows that 7]^^ —7- \ipt){ipt\ in trace-norm, as N ^ 00. 

Remarks. The existence of T > such that p.lUp is satisfied is a consequence of the local 
well-posedness of the nonlinear Hartree equation (jl.Sp . see [IG]. Similar methods to the ones used to 
prove (jl.lip can be employed to show the convergence of higher order marginals 7]^^ with the same 
rate {N~^''^ + aj\f) for any fixed fc G N. If we are satisfied with a slower rate for higher marginals, 
a simple argument, outlined in Section 2 of [19j, shows that (11. lip immediately implies that for any 

The first ingredient of the proof of Theorem 11.11 is the observation that the bound (jl.lOp on the 
ii^^'^-norm of the solution of the Hartree equation (jl.5p . together with the assumption if G H'^(M.'^) 
on the initial data, implies an upper bound on the ii^^'^-norm of the solution of the regularized 
Hartree equation 

^^,4'^^ = x/r^A^l°) - a(^— * U-^\'U-^ , (1.12) 

V I • I + a / 

uniform in the cutoff a > (actually, we prove that \\ipt — (f" \\jji/2 is small, of the order a^'^). 
By the propagation of regularity for the solution of the Hartree equation (both the original equation 

(jl.Sp and the regularized equation (I1.12P ). we also obtain a bound for the norm \\ip). \\h^ uniform in 
a > and in t G [—T,T]. 

The second ingredient in the proof of Theorem 11.11 is the method developed in [21] to establish 
the convergence to the Hartree dynamics for a system of non-relativistic bosons. This method is 
based on the use of a Fock space representation of the many boson system, and on the study of 
the dynamics of coherent states. When analyzing the time-evolution of an initial coherent state, it 



is possible to isolate the Hartree component of the evolution. Moreover, as first observed in [T5] 
(and later in [l3]), the evolution of the fluctuations around the Hartree dynamics can be expressed 
through a two-parameter group of unitary evolutions with an explicit time-dependent generator. The 
problem then reduces to deriving a bound for the growth of the number of particle operator (which 
measures the "number" of fluctuations, after second quantization) with respect to this evolution. 
The crucial observation is that the bound derived in [21] for non-relativistic particles can be easily 



extended to the relativistic setting, once a uniform bound on Hv^J" ||//i is available. 

In Section [21 we show the necessary bounds on the solution ip^' of the regularized equation 
(J1.12p . In Section [3] we introduce the Fock space representation, we define the coherent states, and 
we discuss some of their main properties. Then, in Section ^ we show that the evolution of initial 
coherent states can be approximated by the Hartree dynamics, and we use this fact to conclude the 
proof of Theorem 11.11 

The next theorem is the second main result of this paper. It establishes the convergence of the 
one-particle marginal density associated with the solution of the regularized A/^-particle Schrodinger 
equation to the orthogonal projection onto the solution of the Hartree equation (jl.Sp . in the energy 
norm. The result holds for all t G [— T, T] under the assumption that the H^'"^ norm of the solution ipt 
of (II. 5p remains bounded in [— T, T] (in other words, under the assumption that there is no blowup, 
up to time T). 

Theorem 1.2. Fix ip G i?^(M^) with \p\ = 1 and set iJj^ = p® . Consider an arbitrary sequence 
On > with ajsf — )• and such that N^aj\[ — )• oo as N ^ oo, for some /? > 0. Let ipN,t = e~'^^N^ipN' 
be the evolution of the initial wave function ipj\j- generated by the Hamiltonian hl.^) and let ^\q\ be 
the one-particle reduced density associated with ipN,t- 

Denote by (ft the solution of the nonlinear Hartree equation ( fj.5|) with initial data ipt=o = V'- Fix 

T > such that 

K := sup ll^/Jf ||j:/i/2 < oo . (1-13) 

\t\<T 

Then there exists a constant C = C{k,T, \\ip\\ff2,/3) such that 



Tr 



(1 - A)'/' (tS - \n){n\) (1 - A)^/^ I < C (^ + af) (1.14) 



for allt &M. with \t\ < T , and for all N sufficiently large. In particular, it follows that Xi^ t — ^ Wt){'^t\ 
in energy-norm, as N ^ oo. 

Remark. We believe that the same arguments used to show (J1.14p can be extended to prove the 
convergence (in energy norm) of the higher marginal 7]\r j- To keep the paper readable, we do not 
follow this direction here. Note that a simple argument, similar to the one presented in Section 2 of 
[19j (and mentioned in the remark after Theorem II. ip . shows that (J1.14p implies 



Tr 

for all k gN. 



(1-A.JV4(^(J) |^^^^^^|«A(l_^^jl/4 



^ C^lkiJ^^,+aT 



The proof of Theorem 11.21 is based again on a Fock space representation of the many body 
system, and on the use of coherent states as initial data. As in the proof of Theorem II. H the 



Hartree component of the evolution of the initial coherent states can be isolated, and the dynamics 
of the fluctuations can be written through a two-parameter group of unitary transformations with 
an explicit generator. To obtain convergence in the energy norm, however, instead of controlling 
the growth of the number of particle operator, we need to control the growth of the kinetic energy 
operator with respect to the fluctuation dynamics. Technically, this step (contained in Proposition 
15. 2p is the most challenging part of our paper. In a sense, the fact that we can control the growth of 
the kinetic energy of the fluctuations implies that, although on the N particle level we are considering 
a supercritical regime, after subtracting the (supercritical) Hartree dynamics, the system on the level 
of the fluctuations is subcritical. In Section [5] we prove the convergence to the Hartree dynamics for 
initial coherent states and we complete the proof of Theorem 11.21 assuming Proposition 15.21 to hold 
true. In Section [6l we prove Proposition 15.21 

Theorems 11.11 and 11.21 show that, as long as a bound on the H^'^-norm. of the solution cpt of the 
Hartree equation ()1.5p is available, the evolution of the marginal densities can still be approximated 
by (fit- Next, we ask what happens if the solution ft of (jl.Sp exhibits blowup. Under the assumption 
that ift blows up as t —7- T, for some < T < oo, we show that also the solution of the regularized 
A^-particle Schrodinger equation ()1.9p collapses, if t — t- T, and, simultaneously, A^ — t- oo. The A^- 
particle wave function ^7v,t collapses in the sense that the kinetic energy per particle, which remains 
finite, uniformly in A^, up to time T, diverges to infinity as t — )• T if simultaneously, A^ — )• oo. In 
order to make sure that the solution of the A^ particle Schrodinger equation remains close to the 
solution of the Hartree equation as t approaches the nonlinear blowup time, we have to assume that 
A^ diverges to infinity sufficiently fast. Physically, this condition imposes restrictions to the range of 
many body systems for which the Hartree approximation is valid close to the blowup time. From a 
different point of view (if we think of the number of particles A^ as fixed), it tells us how close to the 
nonlinear blowup time we can expect the Hartree dynamics to be a good approximation for the real 
many body quantum evolution. 

Corollary 1.3. Fix (/? G H'^[M.'^) with \\ip\\ = 1 and set ip^ = f® ■ Consider an arbitrary sequence 
aj\f > with UN — >• and N^a^ —)• oo as A^ —)■ oo, for some /3 > 0. Let V'iv,t = e~^^N^'il:N' be 
the evolution of the initial wave function tp^ generated by the Hamiltonian hl.^) and let 7)^^ be the 
one-particle reduced density associated with ipN,t- 

Denote by (ft the solution of the nonlinear Hartree equation lil.5\) with initial data ipt=Q = f- 
Suppose that Tc > is the first time of blow-up for iff In other words, assume that 

Kt := sup ||(/9s||j:/i/2 < 00 
0<s<t 

for all t < Tc, and 

W'PiWh^/'^ — )• 00 as t — 7- T~ . 

Then, for any fixed t G [0, Tc) there exists a constant Ct > such that 

11(1 - A,jV4^^^j2^(i _ A)l/2^(;) < Ct 

uniformly in N £ 'N. Moreover, for t £ [0,Tc), there exists N{t) £ N with N{t) — )• 00 as t ^ T~ , 
and such that 

||(l-A.jV4^^^^^^j|2^Tr(l-A)V2^(;(\)^^^oo as t ^ T" . (1.15) 

In other words, the kinetic energy per particle is uniformly bounded in N , ifO<t<Tc but it diverges 
in the limit t — )■ T~ , if at the same time, the number of particles tends to infinity sufficiently fast. 



Remark 1. The existence of blow-up for solutions of the nonlinear Hartree equation (jl.Sp has 
been proven in under the assumption that the initial data 99 is spherically symmetric and that it 
has negative energy i?Hartree(</') < 0; see (|1.2|) (this is possible if A > A^;^). 

Remark 2. The fact that 'ipN,t collapses at some point in the interval [0, Tc\ follows already from 
the blow-up of ipt at time T^. and from the Theorem II. 1[ This fact, which was pointed out to us by R. 
Seiringer, follows from the general observation that the kinetic energy of an L^-limit is always smaller 
than the limit of the kinetic energy. This argument, however, does not prove that the collapse takes 
place at time Tc nor that the blow-up of the Hartree equation accurately describes it. 

Proof of Corollary FTTgJ Set f^N = N~'^/'^ + un] then /3n —^ as N -^ 00. For every < t < Tc there 
exists, by Theorem 1.2, a constant Ct, depending on l3,K,t,t, \\(p\\jj2, such that 



Tr 



(1-A)^/V7« 



|^,)(^,|)(l-A)i/^ 



< QA 



TV 



for all < s < t. In particular this implies that 



\l/2^,(l) 



TY(l-Ar/^7]v;,-||v.i||^v. 



<Ct^N 



For < t < Tc, choose now N{t) sufficiently large , so that 7Ar(i) < {Tc — t)/Ct (this is certainly 
possible because Pn — )• as A^ — )• 00). Then 



as t — >• Tc. Since ||(^t||j:/i 



/2 ^- 00 as t —7- Tc, this implies that 

Tr(l-A)V2^W^^^^oo. 
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2 Bounds on solutions of nonlinear Hartree equations 

In this section, we study properties of the solution of the nonlinear Hartree equation (jl.Sp . In 
particular, we need to compare the solution of (jl.Sp with the solution of regularized Hartree equations 
(like ()2.ip . with a > 0). To this end, we first need to establish the property of propagation of initial 
regularity, under the assumption of a bound on the i^^'^-norm. 



Proposition 2.1 (Propagation of regularity) . Fix s > 1/2 and a > 0. Let ip G i/*(M'^) with \\ip\ 
Let ipt denote the solution to the nonlinear Hartree equation 



1. 



idtift = Vl - A(^i - A 



-|- a 



*\W )^t 



(2.1) 



with the initial condition (pt=o = ^- Fix T > such that 

K := sup ||c/9i||j:^i/2 < oo . (2-2) 

\t\<T 

Then there exists a constant v = i>{k,T,s, W^Wh") < oo (but independent of a) such that 

sup \\^t\\H'> < V ■ (2.3) 

\t\<T 



Proof. We follow here the proof of [16^ Lemma 3] with some modifications. Let J(v') := ((I " I + 
a)~^ * |(/9p)(/9. Then we claim that 



In fact, \\J{^)\\h^ < II JMIb + ||(-A)^/VM||2 and 



for all (p G if*' 



s/Tn,3\ 



\\Ji^)h 
Moreover, 

||(_A)*/2j 



< 



I , *lv?nv^ 

• + Q / 



< 



*\ip\' 



\vh ^ IIV'llHi/all'^lb 



(2.4) 



(2.5) 



.A)^/2 



+ a 



*\ip\^ if 



< 



< 



A*I(-A)^/VPI 



Iv'lls + 



+ a 



m 



-A)^/V| 



(2.6) 



i-^Y' \v\ kMMn^^ + Il95|l:^i/2||</'lk= 



^ IIV'llHi/allV'lk^ 



Here we used the generalized Leibniz rule (see Lemma 12. 4p in the first inequality. In the second 
inequality, we used the Hardy-Littlewood-Sobolev inequality and the Sobolev inequality \\ip\\3 < 
1 1 V' 1 1/^1/2 to bound the first term, and Kato's inequality 



sup /d,M^<f /d, |V|VV(,) 
\x-y\ 2 J 



vr 



< 17 IIV'lli:^i/2 



(2.7) 



to bound the second term. Finally, in the third inequality, we used again the generalized Leibniz 
rule. This shows (127 



Next, we write (ft as 



^t 



-Wl-At 



(f + iX dse 



-Vl-A (t-s) 



1 



. + OiN 



* \^s\ ^. 



and we obtain, by ([27iD and ([272]) . that 

IIV'illH^ < llV'lk^ + / \\J{fr)\\H^<dT < \W\\h' + tZ^ I \\<Pr\\H=<^T. 

Jo Jo 



The proposition now follows applying Gronwall's inequality to (j2.9p . 



(2.8) 

(2.9) 
D 
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Next, under the assumption that the solution ipt of the Hartree equation (jl.Sp with initial data 
if G if^(R^) stays bounded in H^'^ in the interval [—T,T], we show the vicinity (in the H^'^-noTva) 
of the solution ipi' of the regularized equation (|2.1|) . with a > and small, to c^^. 

Proposition 2.2. Fix Lp G //^(M^) with \\(f\\ = 1 and let tpt denote the solution of the nonlinear 
Hartree equation / fi.5|) . with initial condition ft=o = f- Let T > be such that 



K := sup ||v9t||j|^i/2 < CXD 
|t|<T 



.(") 



For a > 0, let ip\ he the solution to the regularised Hartree equation 



idt^r = '/^^^vt -^ 



[a] 



1 



(")|2V„H 



+ Q 



m I ft 



with initial condition ipt=o = f- 

Then there exists a constant C = C(T,k, \\^\\h'^) < oo, such that 



m-ft 



Hi 



< C a for all \t\ < T and all a > 0. 



(2.10) 



(2.11) 



(2.12) 



Moreover, if we assume additionally that cp € H^ 
D{T, K, W^pWh^) < oo such that 



k*-9^1"^Lv. < Da'/' 



then we can also find a constant D = 

(2.13) 



for all \t\ < T and < q < 1. 



Proof. We start by proving ()2.12p . Since ^ G H'{M.'^), we can find, by ()2.10p and Proposition 12. H 
ly = i'{T,K, WifWjji) < oo such that 

sup \\(pt\\m < t^ ■ (2-14) 

\t\<T 

Let te [-T,T]. We have 



Therefore 



d II (a)||2 

— \\iPf — ip; 9 

dt"^ V-t 112 



2ARe(^„^(")) 
1 



2Xlm{ipt, [-r-r* \^t\ - , , , 

\x\ \x\ + a 



1 I (")|2\ (a) 
*m I ]ft 






dt 



I (")||2 



< 2|A| 



+ 



if 



a 



*IV5tP)(<^^'' -ft)) 



x\{\x\ + a)) 



(2.15) 



11 



The first summand in tlie r.h.s. of (j2.15p can be estimated as 



{"ft 



a 



\x\{\x\ + a)) 



Wt?){^t^-n)) 



< 

< a 
< 



a 



\x\{\x\ + a)) 



m\ 



m-ft h 



1 



* ml 



m-ft h 



II (")ll II Il2 ^ 2|| (a)|| 



The second summand can be estimated by 



< 



< 



\x\ + a 

dxdy\ipt{x)\ — — \ipt{y) - ip't'\y)\ {\My)\ + I </'!"'' (y) I) I'p't'Hx) - Mx)\ 

;xR3 \x-y\+a 



dxdy|(/9t(x)p 



1 



\x -y\ 



< 



+ 2 



1 



\My) - 4''\y)\' 



dxdy {\My)\' + \^riy)\') \v>ri^) - Mx)\' 



^ ^11 I|2 , i\|| {")||2 ^ n , 2\ \\ {")||2 

where on the last line we used Hardy's inequality 



sup /dyMy^<4/dy|Vv9(y)|2<4||v9||^, 



(2.16) 



Thus, ([2T5]1 gives 



a II {")||2 

^ll^t-V^. II2 



^ /■-t I 2\ / II (a) II I II {' 

< {1 + u )\^a \\ipt - 'Pi 'h + \m - fl 

By Gronwall's inequality, we find C = C{v, T) with 



{")||2 
2 



(")ll ^ r^ 
Pt - 'ft h < <-^" 



(2.17) 



(2.18) 



for all a > 0. 



Next, we prove (|2.13p . To this end, it is enough to show that there exists D = D{T^k,,\\(p\\h2) 
such that 

||(-A)V4(^,_^H)||^ < Da^'\ (2.19) 

Note that, since if £ H'^{W^), we can find, by (j2.10p and Proposition 12.11 u = i>{T,k, \\(p\\h^) < 00 
such that 

sup llv^tllz/a < z^ • (2.20) 



.(") 



\t\<T 

We write ipt and ifl"'' using their Duhamel expansions 

'Pt = e~ 

Jo 
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Jo ^\x\ ' 



"Ps 



and 



(a) 



+ Q 



Jo M^M 

respectively. Thus 

||(-A)V4(^,_^("))||^ < \X\ fds\ (-A)V4(J_*|^^|2)(^^_^H) 

Jo I ^Fl ^ 

.A)V4 



*|J°)|2V» 



'/'s 



+ 
+ 



(-A) 



1/4 



|x|(|x| + Q 
1 



H^.PWi") 



(2.21) 



(Iv^sp-|^(-)P)W(-) 



|x| + a 
Further decomposing the second and third term in the parenthesis we find 

-A) 



+ 
+ 
+ 
+ 



(-A) 



1/4/ 


a 


[ 


\x\{\x\ + a) 


1/4/ 


a 


I 


\x\{\x\ + a) 


1/4/ 


|x|+a ^'^ 


1/4/ 


^ *(h 



* \^s\ y^. 



Iv^i"¥))v^. 



V X + a / 



(2.22) 



The first term is bounded by 

(-A)V4(J_*|^^|2)(^^_^W) 

\\x\ / 



< 



-A)V4/i 



+ 



FT * Wi 
\x\ 



* Ws 



IVs-vi^'^h 



a/4 



i-AY^V^s-^, 



(")) 



(2.23) 



where we used the generahzed Leibniz rule (see Lemma 12. 4p . Next we observe that, by Kato's 
inequality (\2.7h and by (I2.20p . we have |||.|~^ * Iv'sPHoo ^ ll¥'s||^i/2 ^ ^^- This, combined with the 
bound 



\l/4 



1 



x\ 



fsWs ;:; Ws\\hi/2 ^z^" (2.24) 

implies (using also (j2.12p ) that 

(-A)i/V^*|^.|')(^s-V=i")) <Avs-A''^\\hs <i^2„ + z^2||(-A)i/4(^^_^W)||2 (2.25) 

\\x\ / 2 



To prove (|2.24p . we rewrite |.|~^ * \(ps\ = -47r(-A) ^|(/3s| • Then 



(-A) 



1 



1/4/^,1 |2 

x\ 



< (-A)-'-i(^, 



-3/4 1 



le = 11^3/2* I '/'si lie- 
Here Gs, s G (0, 3), is the kernel of the operator (— A)^*' ^ which is explicitly given by 



(2.26) 



(2.27) 
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with C3/2 = 7r^\/2/r(|). From (j2.26p . we conclude by the Littlewood-Hardy-Sobolev inequahty that 



(-A)^/^(^ 



Ws 



|2 



< IIG3/2 * Iv'sPlle < \Ws 

and thus (l2:2iD follows. 

The second term on the r.h.s. of (j2.22p is estimated again by the generalized Leibniz rule as 



< 



(-A)V4(. 



a 



rr ( x + a 



* IV^s 



Iv'.-v'i'^lb 



(2.28) 



+ 



a 



Since 



-A) 



1/4/ « ^ 1,^ |2 

|x|(|x| + a) 



XX +a, 



< 



* \^s 



a 



XX + q) 



|(-A)V4(^^_^H)||^. 



\{-^Y"Ws?\ 



< aV2||(_A)VV||3||^.||6 < a'^' W^^Wl. < a'/' u' 



(2.29) 



and, by (f2l6]) . 



a 



|x|(|x| + a) 
we find, using (j2.12p . that 



Ws 



< a 



T2*\'P^\ 



^ "llv'sll/fi ^"'^^ 



(-A) 



1/4 



a 



X ( X + a) 



(2.30) 



*|^.P)(v^.-c^i")) ^ < a3/V + a^l(-A)V4(^^_^H)||^. (2.31) 



The third summand in (j2.22p is estimated (again using (|2.16p ) as 



-A) 



1/4 



a 



x|(|x| + a) 



* \Vs\ •^. 



< 



a 



X ( X + a 



+ 



-A)V4^ 



* \Vs 



a 



\x\i\x\ +Q, 



* W^ 



A)i/Vll2 



where we used ()2.30p . the Sobolev inequality Hv'sHe ^ l|y's||_H'i ^-nd the bound 



(-A) 



1/4. 



a 



X ( X + a) 



Ws 



< a J-*|(-A)i/Vsn 



a||(-A)i/^((^V=,)||3/2 

" II'/'s||hi||v5s||//i/2 



< 



< 



(2.32) 



where we used the Littlewood-Hardy-Sobolev inequality, and, in the last inequality, the generalized 
Leibniz rule (see Lemma 12.41 below). 
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The fourth summand on the r.h.s. of (j2.22p is bounded by 



(-A)^/^fo^ 



\x\ + a 



*(|^,|2-|^(")|2m^. 



< 



(_A)i/4 1 



x\ + a 



+ 



1 



-{\^s\'-\^t\') 



*{Ws?-WT^V) J|(-A)vv.||3 



2+e 
\l/4, 



\^s\\ -2{2+e) 



(2.33) 



|x| + a 
for arbitrary e > 0. The second term on the r.h.s. of (j2.33p is estimated by 



^*(|V..P-Iv^i°)p) J|(-A)VV.||3 < 119^.^1 



\x\ + a 



\x\ ' ' 



6/5 (2.34) 

(")|lo 



where we used the Sobolev inequahty on the first, the Hardy-Littlewood-Sobolev inequaUty on the 
second, the Holder inequahty on the third, and, finahy, again the Sobolev inequality in the fourth 
line. 



As for the first term on the r.h.s. of ()2.33p . we notice that 
1 



(-A) 



1/4. 



\x\ + a 



* [I'Ps 



lv^i"¥) 



2+e 



l'Psl| 2(2+e) 



< 



< 



(-A) 
(-A) 



1/4. 



1 



1/4 



\x\ + a 

1 
\x\ + a 



2+e 



2+e 



^,-(^(-)|(|(^,| + |v9(-)|)||J|^, 



Wi 



s\\m 



(2.35) 



^,-y,<f)\\,{\\^,\\, + \\^f)\\,) \\^ 



(")| 



s\\m 



< i^a'-' 



where in the last step we used the bound 



(-A) 



1/4. 



1 



\x\ + a 



a 



2+e 



for all e > 0. The bound (12.360 follows from the pointwise estimate 

1 



f(_A)V4^_') 
V \x\+aj 



< 



valid for all x G M^. To show (j2.37p . we observe that 



(-A) 



1 



{\x\ +a] 



(x) 



and therefore 



-A)i/S 



1 



\x\ + a 



(-A) 



-3/4 _ 



a 



\x\(\x\ + a) 



(|x|+a)3/2 



2a 

|x|(|x| + a) 



< dy 



(2.36) 



(2.37) 



a 



y|3/2 \y\{a + \y\y 



(2.38) 
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We assume first that \x\ > a. From (j2.38p we find 



(-A) 



1/4. 



1 



\x\ + a 



< 



< 



dy 



a 



+ 



\x-y\>\x\/2 \x - 2/P/2 \y\ia + |y|)3 y|^._y|<|^|/2 \x - y|3/2 \y\{a + \y\)^ 



dy 



a 



dy- 



a 



+ 



dy- 



|x|3/2 7 |y|(Q, + |y|)3 |x|3 7|^,_j^|<|^|/2 |x-y|3/2 

where we used the fact that |x — y| < \x\/2 imphes \y\ > \x\/2. Expficit computation hnphes that 



-A) 



1/4. 



1 



\x\ + a^ 
For \x\ < a we notice that, by (|2.38p . 



(x) 



< 



< 



|x|3/2 - (|x|+a)3/2 



-A) 



1/4. 



1 



\x\ + a 



dy 



a 



\x~y\>a |x-y|3/2 |y|(a + |y|)3 



+ 



i 



for all \x\ > a. 



dy 



(2.39) 



a 



\x-y\<a |x-2/|3/2 \y\(a + \y\y 



<^ I dy , ,. . , „, + — 



a 



3/2 



dy' 



y|(a + |y|)3 a^ 7|3._y|<„ \x - y\^/^\y 



(2.40) 



Since |x| < a and |x — y| < a imply that |y| < 2a, the last term is bounded, for |x| < a, by 



/ 

J\x- 



dy- 



< 



dy 



+ 



1 



'\x-y\<a \X - yP/2|y| ^ J\x-v\<a \x - yp/^ J\y\<2a l^p/^ 

Inserting back in (j2.40p . it follows that 



<r.l/2 



(-A) 



1/4. 



1 



\x\ + a 



[x] 



<J_< I 

~ a3/2 ~ (|x|+a)3/2 



for all \x\ < a. 



Together with ([239]) . this implies (fOT]) and therefore (fOSj) . Combining (j23i]) with (fOSj) . we 
obtain the bound 



(-A)V4(^_ , (|^^|2 _ |^W|2))^^ < z.||(-A)V4(^^ _ ^H)||2 + ^ai- 

for all e > 0. 

The fifth summand in ()2.22p is estimated as 

(-A)^/^(^ * (I^.P - |9^i")p))(^. - ^i")) , 
\|a;| + a / 2 

((-A)V^^)*(|,.P-|.i"f) 



(2.41) 



< 



l¥^s-¥^i°)||2 



+ 



1 



X + a 



* {\y^s\' - l^'fY) ||(-A)V4(^^_^(")) 



(2.42) 



< 



-A) 



1/4. 



1 



Ws - ^, 



\x\ + a 
< aV2 + ||(_A)i/4(^^_^W) 



i")|li + -ll'^.-'^i"^ll2||(-A)V4(^.-^W) 



where we used the generalized Leibniz rule in the first inequality and the bounds ()2.12p and ()2.37p 
in the last inequality. 
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Inserting the estimates (l2:25|) . (ITSTT) . (f2:32|) . ([MI]), and (i2:i2]) into (i2:22]l yields (using that 
z/ > 1 and the assumption a < 1) 



|(-A)V4(^,_^W)||^ <-^^'d.|||(-A)V4(^^_^W)||^ + «i/2|. 



(2.43) 



Eq. (j2.13|) fohows by Gronwah's lemma. D 

In the next corollary, we summarize the consequences of the bound on the iJ^' ^-norm of ift (and 
of the assumption y? G H'^(W^) on the initial data), that will play a crucial role in the many body 
analysis. 

Corollary 2.3. Fix s > 2 and tp G H'^i^). Let ipt and, for any q > 0, 93^° be the solutions of the 

nonlinear Hartree equations hl.5\) and, respectively, \2.11\) with initial data ^pt=Q = vt—Q — V (y^t is 

the maximal local solution of lil.5\) in i7^'^(]R^); (f^ , on the other hand, is known to exist globally 
in H^/'^{R^)). FixT >0 such that 

K := sup ||(/3i||j:^i/2 < 00 . 
\t\<T 

Then there exists v = v[s,T,k, W'^Wh") < 00 independent of a such that 

II (")ii ^ 
sup Wfl \\h= < z^ 

\t\<T 

for all a > small enough. 

Proof. Since s > 2, Proposition 12.21 implies that 

II (") II ^ n 

sup \\(pl 11/^1/2 < 2k 
|t|<T 

for sufficiently small a > 0. The claim follows then by Proposition 12. li D 

To conclude this section, we state the generalized Leibniz rule for fractional derivatives. For a 
proof of this lemma, see 



Lemma 2.4 (Generalized Leibniz Rule). Suppose that 1 < p < 00, s > 0, a > 0, /3 > 0, and 

1/pi + 1/qi = 1/p with i = 1,2, 1 < qi < 00, 1 < Pi < 00. Then there exists a constant c = 
c{p,pi,p2,s,a, f3) < 00 such that 

U-AY/'{fg)\\, < c(||(-A)(«+°)/2/lkll(-A)-"/25||,, + \\{-A)-^/'f\Ui-A)(^+(^y'g\Q 

for all measurable functions f,g for which the r.h.s. is finite. 
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3 Fock space representation 

In this section, we introduce a Fock-space representation of our system, and we define coherent states. 
The bosonic Fock space over L^(]R^, dx) is defined by 



jr = l2(M^ dx)®»" = C e L^(]R3n^ dxi . . . dx„) , 

n>0 ra>l 

with the convention L^(]R^)®^'^ = C Vectors in J^ are sequences t/j = {il^^"''}n>o of n-particle wave 
functions V' £ L^(M^"). On J", we introduce the scalar product 

It is simple to check that, with this inner product, J-" is a Hilbert space. States with N particles and 
with wave function ipN G Lg(M^ ) are described on T by the sequence {V' }n>o where V' = for 
all n 7^ A^ and ip^^' = V'at. The vector {1, 0, 0, . . . } G J-" is called the vacuum, and will be denoted 
by 17. 

The number of particles operator M acts on J^ according to {Mil))^"^' = nip^^^ for all n G N. 
Eigenvectors of J\f are vectors of the form {0, . . . , 0, t/:^"^' , 0, . . . } with a fixed number of particles. 

For arbitrary / G L^(M'^) we define the creation operator a*(/) and the annihilation operator 
a{f) on T by 

1 " 

(a*(/)^)^"'' {xi,---,Xn) = ^Y] f{Xj)i)^'^~^'>{xi, ..., Xj-i,Xj+i, ...,Xn) 

^ i=i (3.1) 

(a(/)^)(")(xi,...,x„) = ^/^m/dx7MV^"+i)(x,xi,...,x„). 



The operators a*{f) and a{f) are unbounded, densely defined, closed operators. The creation oper- 
ator a*{f) is the adjoint of the annihilation operator a(/) (note that by definition a{f) is anti-linear 
in /), and they satisfy the canonical commutation relations 

Hf),a*{g)] = {i\g)L^i^,), [a{f),a{g)] = [a* {f) , a* {g)] =0. (3.2) 

We will also make use of operator valued distributions a* and Qx (x G M^), defined so that 



«*(/) = / dx/(x)a* 

(3.3) 



«(/) = / dx/(x)aa, 

for every / G -L^(M^). The canonical commutation relations assume the form 

[a^, a*y] = 5{x - y) [a^, ay] = [a* , a*] = . 
The number of particle operator, expressed through the distributions axja^, is given, formally, 

by 

M = dx a*i.ax ■ 



The following standard lemma provides some useful bounds to control creation and annihilation 
operators in terms of the number of particle operator N . 
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Lemma 3.1. Let f e L^(M'^). Then we have 

HfM < 11/11 ||AAi/2^|| and KUM < 11/11 II (M + 1)1/2 _^|| ^3^4) 

For an arbitrary ip ^ J-, we define the one-particle density 7^ associated with ip as the positive 
trace class operator on L2(M'^) with kernel given by 

^i'^ (^; y) = i^^^j^^) (^' "y^-V') • (3.5) 

By definition, 7^ is a positive trace class operator on L2(M^) with Tr 7! = 1. For every A^-particle 
state with wave function ^tv G Lg(R^^) (described on T by the sequence {0, 0, ... , ipN,0, 0, . . . }) it 
is simple to see that this definition is equivalent to the standard definition. 

For any sequence a = {a^), with a^ — )• as A^ — )• 00, we define the Hamiltonian Ti'^ on J^ by 

(7^°^)W = (-^^)(n)^(«), with 

Using the distributions ax,a'^, 11% can be rewritten, formally, as 

Ti^ = / dx a* (1 - A^)i/2 Ox - 77T7 / dxdy -^ t- 0*0*0^0^, . (3.6) 

J 2N J \x - y\ + ON '^ 

By definition, the Hamiltonian 'H'^ leaves sectors of J-" with a fixed number of particles invariant. 
Moreover, it is clear that on the A^-particle sector, 'H'^ agrees with the Hamiltonian H^. We will 
study the dynamics generated by the operator 'H'^. In particular we will consider the time evolution 
of coherent states, which we introduce next. 

For / G L2(R^), we define the Weyl-operator 



T^(/) = exp(a*(/)-a(/)) (3.7) 

and the coherent state V'(/) S T with one-particle wave function / by V'(/) = W{f)Q,. Notice that 



"' ^n vn! 

n>0 n>0 

where f^"' indicates the Fock-vector {0, . . . , 0, /®", 0, . . . }. This follows from 

exp(a*(/) - a(/)) = e-H^H^^ exp(a*(/)) exp(-a(/)) 

which is a consequence of the fact that the commutator [a(/),a*(/)] = ||/|p commutes with a{f) 
and a*{f). From Eq. p.Sp we see that coherent states are superpositions of states with different 
number of particles (the probability of having n particles in ^p{f) is given by e""-^" ||/|p"/n!). 

In the following standard lemma we collect some important and well known properties of Weyl 
operators and coherent states. 



Lemma 3.2. Let f,g ^ L 
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i) The Weyl operators satisfy the relations 

Wif)Wig) = VF((7)Ty(/)e-2^^™(/'9) = W{f + g)e-'^"'^f'3) . 

a) W{f) is a unitary operator and 

wifr = wif)-' = w{-f). 

Hi) We have 

W*{f)a,W{f) = a, + f{x), and W*{f)alW{f) = a*, +J{x) . 

iv) From Hi) we see that coherent states are eigenvectors of annihilation operators 

a^iPif) = f{x)^{f) ^ a{g)^{f) = {g, f) L^^^if) . 

v) The expectation of the number of particles in the coherent state ^(/) is given by ||/|p, that is 

{il^{f)Mi'{f)) -- 



Also the variance of the number of particles in ipif) is given by ||/|p (the distribution of M is 
Poisson), that is 

{^{f),M^^{f))-{i^{f),M^{f)?-- 



vi) Coherent states are normalized but not orthogonal to each other. In fact 

{^{f)M9)) = e-Kll/P+ll9lP-2(/.s)) ^ KV(/),V'(5))I = e^^ll^-^ll' . 

4 Time evolution of coherent states and proof of Theorem 11.11 



In this section we study the time evolution of an initial coherent state ijj{y Nf) = W{yNf)VL, for 
</? E H'^i^) with ||c^|| = 1. The expected number of particles in the coherent state ip{\/Nip) is A'". 
Therefore, we may expect the evolution generated by T-L'^ on ^(viVt/?) to have a mean- field character. 
In particular we may expect that e~ 'Nip[\/N(p) ~ il){\/N(ft) where ipt solves the nonlinear Hartree 
equation (ll.Sp . We will prove that this is indeed the case, under the assumption that ipt remains 
bounded in /7^'^(R^) in the time interval [—T,T]. 



Theorem 4.1. Fix ip G i?^(M'^) with ||(/9|| = 1 and an arbitrary sequence a^ > such that ajsf — )• 
as A^ — )• oo. Let ip{N, t) = e~''^'^NW {\/lVip)i} be the evolution of the initial coherent state W{VN(f)fl 
generated by the Hamiltonian \3.0) . Denote by Tjy ^ the one-particle reduced density associated with 

i;{N,t). 

Let ift be the solution of the nonlinear Hartree equation iL5\) . with initial data ipt=o = ^- Fix 
T > so that 

K := sup llv'tll/fi/a < oo . (4.1) 

\t\<T 

Then there exists C = C{T,k, \\<f\\H2) < oo such that 



Tr r^l- \'Pt){Vt 
for allteR with \t\ < T. 
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<C{^ + aN 



Proof. The proof of Theorem 14. II is analogous to the proof of Theorem 3.1 in |21j . For completeness 
(and because some of these arguments will be used later on), we explain here the main steps. 
Since \ipt){^t\ is a rank one projection, it is enough to show that 



^(1) 

N,t 



\ft){<^t 



< C I — + OAT 
HS V A^ 



where ||^||hs ~ TvA*A is the Hilbert-Schmidt norm of A. This follows from the remarHll that the 
operator Fj^ ^ — \ipt){ipt\ can only have one negative eigenvalue. Sine the trace vanishes, the absolute 
value of the negative eigenvalue must be the same as the sum of all positive eigenvalues. For this 
reason, the trace norm is twice the operator norm, which is of course bounded by the Hilbert-Schmid 
norm. 



Suppose now that (p^^ denote the solution of the regularized Hartree equation 



idtifl 



(oiv) 



^r^Avpj' 



("iv) 



1 



. + aN 



Hv^^^pW^) 



(4.2) 



with initial data v^J^g ~ ^- ^^ Proposition 12.21 (see, in particular, ()2.12p ). and since 



m )m \-m){^t 



HS 



< 2\\ipt -Lfl 



(a]v)| 



it is enough to prove that 



\^ N,t Wt /\Vt I 



|HS 



c 

< — 

- N 



(4.3) 



for a constant C depending on T, k, ||(/9||ji^2. In order to prove (j4.3p we write the difference of the 
kernels of F^yj and \ipi''^ ){ip\"^ \ as (compare with (3.4) in [21j) 






N 



N 



n,W*iVNip)e''^^\a^ - VNipf'"{x))e-''^^'W{VN^)n 



("iv)/ 



iV-Ni^ 



N 



^,W*{VN^)e^^^\al - VNlp^^''\y))e''^^*W{VN<p)Vt 



]- (vL,W*{VNip)e^^^\a*y - ViV^;"^)(y))(a^ - VN^^^''\x))e-'^^'W{^ip)Vt^ . 

(4.4) 



Then, following (3.5)-(3.8) in [21], we can define the two-parameter group of unitary transformations 
UN{t; s) by the Schrodinger equation 



idtUjyit; s) = CN{t)hlN{t', s) and Un{s; s) = 1 



(4.5) 



^We learned this argument from R. Seiringer 
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with the generator 



CNit)= [dxalil-A,)'/^a.-X [ dx ( --^ * Iv'^'^^rt (^) «>- 

J J \\.\+aN J 

-A [dxdy. ^ (^("-)(x)(^^-)(y)aX. 



J \x - y\+ aN ' ^ ^ " " "" ^ 

- ^ / dxdy \-— Ur^ (x)(^J"-) (yXa; + c^^"-) {x)lpl'^"-^ {y)a,ay) (4.6) 

^ J F - 2/1 + Oat V ^ / 

VN J \x-y\+aN ^ " > 

- T^j \ d2;dy — ala*aya^ . 

2N J \x — y\ + UN ^ 

It was observed by Hepp in [15] and then by Ginibre-Velo in [13] that 

UUt;^)axUN{t]^) = t^*(\/iV(/.)e^^^*(a^ - \/iV(^J"^)(x))e-^^^*VF(\/iV(/p) . (4.7) 

Therefore it follows from (j4.4p that 



T%{x,y) - ^t^\x)-^t''\y) = T7 {nMN{m*ala^UN{t;m) 



1 

iV 

^^-)(x) 



+ ' r^ ' {n,UN{t;Ora;UN{t;0)n) (4.J 

+ ^* ^'^^ {n,UNit;0)*a.MN{t;0)n) . 



N 

To get an optimal bound on the error, we also introduce, similarly to (3.9) and (3.10) in ^21j, the 
modified evolution U]\[{t; s) defined by the equation 

idtUN{t;s) = CN{t)UN{t;s) with Un{s]s) = 1 (4.9) 

with the time-dependent generator 

CN{t)= [dxal{l-A,f/^a.-X I dx (--^ * Iv't"'^^! (^) <«x 

J J \\-\+OLM J 



X [dxdy- l^^[^-\x)^l^-\y)a;a, 

J \x-y\+aN ^ 



A /• , , 1 



dxdy ^ '—— U"-)(x)^^-)(y)a:a; +^f-)(x)c^f-)(y)a.aj 

\x — V\ + (Xn V 



2 y \x — y\+ On 

~ ^7^ / da^dy -. r- ala*aya^ . 

2N J \x -y\+aN " 

Since Un commutes with the parity operator (—1)-^, we have 

'n,UN{t;0)*ayUNit;o)n) = (n,UN{t;0)* aiUN{t; o)n) = 0. 



(4.10) 
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Therefore, we can write 
= ^{n,UNit;OyalaMN{t;0)n) 



+ 



TV 



n,u*N{t; o)a^ (uN{t; 0) - i(N{t; 0)) n) + (n, {u*j^{t- o) - w^{t- 0)) ajANit; m 



which leads, after multiplying with a Hilbert-Schmidt observable J and taking the trace, to the 
bound 



TYJ r 



(rS-l^^^)(^^^l) 



< ^^ {UN{t;f))n,MUN{t;f))n) 



N 
^II-'IIhs 



TV 
^Ik IIhs 



\\{UN{t; 0) - UnH; 0))fl|| || {M + lfl^UN{t; 0)O|| 
\\{UN{t; 0) - UN{t; 0))f]|| || {M + l)^'^UN{t; 0)O|| . 



(4.11) 



To conclude the proof of the theorem, we combine the last bound with Proposition 14.21 Proposi- 
tion [4]3l and Proposition 14.41 below. D 

The next proposition shows that expectations of powers of the number of particle operator, 
evolved with respect to the fluctuation dynamics Un-, stay bounded up to time T. Note that to prove 
Theorem 14. H it would be enough to have (|4.12p for k = 1 and for ip = Vl; for later use, however, it 
is useful to consider arbitrary A; G N and ip £ T). 



Proposition 4.2. Suppose that the assumptions of Theorem \4-l\ are satisfied. Suppose moreover 
that the unitary evolution Ui^{t\s) is defined as in |^.5[ j. Then, for every /c G N, there exists 
C = C{k,T,K,\\(p\\H2) such that 



{UNit;0)^,M''UNit;0)^) <C{ij,{M + l 



,2fc+2 



^) 



(4.12) 



for every ip £ T and every t G M with \t\ < T. 



A similar estimate is also needed to control the growth of the expectation of the number of 
particle operator with respect to the modified dynamics Un introduced in (j4.9p . 



Proposition 4.3. Suppose that the assumption of Theorem \4.1\ are satisfied. Suppose moreover that 
the unitary evolution U]\f{t; s) is defined as in |^.g[). Then there exists C = C{T, k, ||95||j|^2) such that 



{UN{t;0)n,M^UN{t;0)n) < C 
for every ip £ T and every t G M with \t\ < T. 



Finally, we need to show that, in the second and in the third term on the r.h.s. of ()4.1ip . it is 
possible to extract one more factor N~^''^ from the difference between the two evolutions. 



23 



Proposition 4.4. Suppose that the assumption of Theorem \4.1\ are satisfied. Suppose moreover that 
the unitary evolutions UN{t; s) andUN{t]s) are defined as in ( [^.Jp and in ^-S^ - Then there exists 
C = C{T,K,\\ip\\^2) such that 



UN{t;0) -UN{t;0)) n 



c 

< 



N 



The proof of these three propositions can be obtained in the exact same way as the proof of 
Proposition 3.3, Lemma 3.8 and Lemma 3.9 in [21]. This fohows by the observation that, on the one 
hand, the kinetic energy (given by the second quantization of the dispersion (1 — A)^'-^), which is 
the only term in the generators C]\f{t) and Cj\i{t) which differs from the generators in [2l] . commutes 
with the number of particle operator (and with all its powers). The other important remark is 
that by the assumptions in Theorem 14.11 (in particular, by ()4.ip ). and by Corollary 12. 3^ there exists 
u = u{k,T, W^Wh^) < oo independent of a^ such that 

II {"iv)ll ^ 

sup Wfi \\m < '^■ 

\t\<T 

The uniform bound on the f/^^-norm of ipt is the only property of cpt that is used in the proof of 
Proposition 3.3, Lemma 3.8 and Lemma 3.9 of [21j. 

Note that the main idea in the proof of Proposition 14. 21 is the introduction of yet another modified 
dynamics WAf(t; s) defined by 

idtWNit; s) = MN{t)WN{t; s) with Wn{s; s) = 1 for all s G M, (4.13) 

with the time-dependent generator 

MNit)= /dxa*(l-A,)V2a^-A [ dx ( —^ * {^'"''^A (x) a^a, 

-X [dxdy-. \- W^"-)(x)<^;"'^)(y)a;a, 

J \x-y\+aN ^ 

- ^ [dxdy- \-— (y,i"-)(x)(^l"-)(y)a:a; + W^"-)(x)^(°-)(y)a..a,) (4.14) 

2 J \x -y\ + aN \ -^ / 

-^ [dxdy- l^a:Ur\y)lM{^f)a;+Tp-t^--\y)aylM{^f))a, 

VN J \x-y\+aN ^ ^ ) 

~^^ ^^^y I n ala*aya^ . 

2N J \x — y\ + UN 

where, for every M > 0, 1m (s) = 1 for s < M, and 1m (s) = if s > M {1m is the characteristic 
function of (— oo, M]). At the end M is chosen as M = const • N. One of the main steps in the proof 
of Proposition 14.21 is a bound for the growth of the expectation of the number of particles w.r.t. the 
cutoffed dynamics WAr(t; s). We state this result explicitly, because similar ideas are used also in the 
next section for the proof of Theorem 15.11 The proof of the next lemma is analogous to the proof of 
Lemma 3.5 in [21]. 

Lemma 4.5. Suppose that the assumptions of Proposition \4^ are satisfied. Let Wn be defined as 
the unitary evolution ^.13 ) with generator Iji4.14^ and with M < const • N. Then, for every k £ N 



there exists C = C{const, k,T,K, ||v^||/f2) such that 

{WN{t;0)n,Af''WN{t;OM < C 
for allt£R with \t\ < T. 
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Proposition 14.21 also allows us to conclude the proof of Theorem 11.11 by writing the factorized 
initial data as linear combinations of coherent states. We follow here the proof of Theorem 1.1 in 
|21j very closely; for this reason, we only discuss the main ideas. 



Proof of TheoremUJi For ipN = V^®^ G L^ 



we write (see Lemma 4.1 in [2T 



VNl Jo 27r 



(4.15) 



with the constant 



d-N 



W 



]\fN/2^-N/2 



Ari/4 



(1) 



The kernel of the one-particle reduced density 7]^ ^ associated with -(/'at 
given by 



{l)r^...^_d% f^^d9i f^^dO. 



N 



i'Nti-;y) = -^ ^ -^e 



2tt 



2 -iBiN ie2N I 



27r 



{W{e-"'WN^)n,a*Jt)a^{t)W{e~'''^^N>f)VL) (4.17) 



-iH'^t 



iB2 



(4.16) 
■^AT is therefore 



where we introduced the notation ax{t) = e^^^^a^e ^"^^^ . As in (4.5)-(4.7) of [2T], we find 
7gi(:r;y) -#-)(,)^h)(x) = f f ^ ^ ^^ ^ ^-^'^^^^'^"^ (^H^; 0)^^, <a.^^^(t;0)f] 

J ,_("iv)/™\ 1-2-K JQ , 

1^ ^{uf,mn,a;^^^^-'^ 



+ 



'N JO 



N 



(4.18) 



where (p®(^~i) actually denotes the vector {0, . . . ,0, (/?®(^~-'^\0, . . . } G F and where W^(t;0) is 
defined as the unitary evolution in ()4.5p . with 97^°^ replaced by e ^p^^ (it is important to observe 
that if (Pt"^ solves the nonlinear Hartree equation, also e (ff^ is a solution). Therefore, we 
conclude that 



Ji) 



^""n) I ^-{oiN) I 



dxdyi7]v;(x;y)-^r^(^)^r^(y) 



< 2 



27r 



d^i 



"TV 

A^2 /p 2tt Jq 27r 



27r 



d^o 



||AAi/X^(t;0)Of ||AAi/X'(t;0)17f + - / dx\fN{x)\ 



N 



(4.19) 



with 



/^(:,) = d% fj ^ (i^^(i; 0)n, a. v^®(^-l) 



(4.20) 



Proceeding exactly as in Lemma 4.2 of ^21j, we find a constant C = C(T, k, W^Wh^) such that 

dx\fN{x)\^<C 



uniformly in N. Proposition 14.21 implies therefore that there exists C = C{T,k, \\<f\\H'^) with 



HS ^/N 



25 



Therefore, it follows from Proposition 12.21 (in particular, (j2.12p ) that 



.(1) 



77v,i - I'/'tXy'tl 



HS 



< c 



+ OAf 



Since |y3t)(y3t| is a rank one projection, and since Tr7]yj = Tr |(/9()((/9(| = 1, the trace norm of the 
difference 7)^^ — |c/?t)(c/?t| is at most two times its Hilbert-Schmidt norm. This completes the proof 



of the theorem. 
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5 Convergence in energy 

We study again the evolution of initial coherent states in the Fock space. This time, we establish 
the convergence of the one-particle reduced density towards the solution of the Hartree equation in 
the energy norm. As a consequence, we obtain a proof of Theorem 11.21 



Theorem 5.1. Fix y? G Il'^i^^ with \\ip\\ = 1 and a sequence on > such that ajv — ^ and 

N^un —7- 00 as N ^ 00, for an appropriate /? > 0. Let ip{N,t) = e~^^'^'NW{\'N(p)Q be the evolution 
of the initial coherent state W{\nVip)Q generated by the Hamiltonian i fg.6]) . Denote by Fj^ ^ the 
one-particle reduced density associated with 'ilj[N,t). 

Let ipt be the solution of the nonlinear Hartree equation M.5\) with initial data ipt=o = ^- Fix 
T >0 so that 



K := sup ||99j||j:^i/2 < 00 . 
\t\<T 

Then there exists C = C{j3,T,k, \\^\\h'^) < 00 such that 

;i-A)V4(rS-|^t)(^t|)(i-A)V4 

for allt£R with \t\ < T. 

Remark. From (j5.2p we can conclude, using arguments similar to the ones used below in the 
proof of Theorem 11.21 (starting from Eq. (IS.lOj) ). that 



HS 



1/2 



(5.1) 



(5.2) 



Tr 



(l-A)V4frg) |^,)(^,|)(1_A)V4 



Ao^n) 



1/2 



Proof. Denote now hy (p). ^ the solution of the regularized Hartree equation (j4.2p with initial data 



^1=0 ~ 'P- Since 



(1 - A)V4 (|^,)(^,| _ |.^("-))(^(°-)|) (1 - A)V4 



HS 



^ II (o'n) II 



and using Proposition [22] (see, in particular, (|2.13p ). it is enough to prove that there exists a constant 
C = C{k,T, \\(p\\h^) such that 



(i-A)V4(ri^)_l^("-))(^("-)|)(i_A)i/4 



HS ^/N 



(5.3) 
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To show (j5.3p , we use again the representation (|4.8p for the kernel of Tj^ ^ , which imphes that 

(1 - A)V4 (rW - |^r-))(^r-)|) (1 - A)V4)(^,y) 

= ^ (f^,Z^jv(t; 0)*(1 - Aj,)V4a* (1 _ A,)i/4 a,ZYiv(i; 0)f^ 



^ (^ A)i/V^ ^ (x) /5^^^^(^.o)*(i _ Ay)'/^a:UNit;0)n 



i'-^Hx) /^,. .. „.... . .,M ... \ (5-4) 

N 



With a Schwarz inequaUty, we find 



fl,ZY^(t;0)*(l - A,)'/^a,Z^7v(t;0)f^ 



HS 



1 I|(P^"^^|P ^^-^^ 

where we defined 

/C = y"dx(l-A,)i/4<.(l-A,.)'/V (5.6) 

to be the kinetic energy operator. The theorem fohows now from Proposition 15.21 below. D 

The key point, in the proof of Theorem lS.H and also in the proof of Theorem 1 1.2 1 is the following 
proposition, which controls the growth of the expectation of the kinetic energy with respect to the 
fluctuation dynamics Un ■ 

Proposition 5.2. Suppose that the assumptions of Theorem 15. Jl are satisfied. Suppose moreover 
that the unitary evolution UNit; s) is defined as in |^.5p . Then there exists C = C{T,k, \\(p\\h'^) such 
that 

{UN{t;0)n,ICUN{t;0)n) < C (5.7) 

forteR with \t\ < T. 

The proof of Proposition 15.21 is given in Section [6j The bound on the growth of the expectation 
of /C can also be used to conclude the proof of Theorem 11.21 

Proof of Theorem \1.2[ Using the representation ()4.18p . we obtain 

((1 - A)V4(^(;) _ |^("-))(c,("-)|)(l - A)y^)ix,y) 

^^ r I "'^^''^^' "^^ " ^-)'^'«-^^'(*; 0)^11 11(1 - A,)i/^a,<^(t;0)f7|| (5.3) 



+ ^ 1(1 - A.)i/vi"-)(x)| rv 11(1 - /\yY'\u'^{t,m\ 



where lA^{t;Q) is defined as the unitary evolution in (|4.5p . with ipf^ replaced by e^^tpf^ . Taking 
the square and integrating over x, y, we find 



dxdy 



(1 - A)V4(^g) _ |^l"-))(c^J"-)|)(l - A)V4)(,,y) 



< 

~ AT 



]-{ [^d9{Uf^it;0)n,ICU^^{t;0)n)] +^ [^de{uf^it;o)n,icu^^{t;o)n) 
^ V JO / vN Jo 



(5.9) 
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Proposition 15.21 implies that there exists C = C{T,k, \\^\\h^) such that 

|(i-A)V4(^(;)_|^(".))^^(".)|)(i_^)i/4|i - c 

By Proposition 12.21 we obtain therefore 



< 



HS ~ A^V4 



Note also that ()5.8p implies that 
dx 



HS 



< C 



1 1/2 

iVl/4 + "^ 



1 / " (aiv) 



< 



JVl/4 I ll'^t 

and therefore, by Proposition 15.21 



2tt 



H1/2+ I de{u]if{t;0)n,icu"N{t;o)n] 





Tr(l-A)V4^^\)(l-A)V4 
Again, Proposition 12.21 implies that 



(l-A)VVi' 



1/4 {oiv)||2 



~ iVl/4 



Tr(l-A)i/SS(l-A) 



a/4 



fi - A)i/4 



w 



<c 



1/2 



]^i7i + "^ 



Last equation, together with (|5.1ip . implies that 

(1 - A)V4^(;) (1 _ A)V4 1(1 _ A)V4^,)((1 - A)V4 



W 



TV[(1 - A)V4^W (1 _ ^)i/4] 11(1 - A)i/4<^,||2 

(1-A)V4(^(;)_|^,)(^,|)(1_A)V4| 



HS 



< 



IHS 



+ 



+ 



Tr[(l-A)V4^S;_)(l-A)V4] 
|Tr[(l - A)V4^a) (1 _ A)V4] _ 11(1 _ A)V4 



v^iiiil 



TY[(l-A)i/47]^,(l-A)i/4] 



(1), 



< C 



1 1/2 

Ari/4 + "^ 



(5.10) 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



(5.15) 



On the I.h.s. of (j5.15p we are now comparing a density matrix with a rank-one projection. The trace 
norm of their difference is at most twice the corresponding Hilbert-Schmidt norm and thus we find 



Tr 



(1 - A)V4^« (1 - A)V4 (1 „ A)V4|^,)(^,|(i _ A)V4 



TV[(1-A)V4^«(1-A)V4] 



(l-A)V4^t||i 



< C 



1 , 1/2 



Ari/4 



(5.16) 
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Combining last equation with (|5.14p . we finally obtain 



Tr 



(l-A)V4(^W_|^^^^^^|)(l_^)l/4 



<||(l-A)i/VfTr 



<||(l-A)VVfTr 



:i - A)V4^« (1 - A)V4 (1 _ A)V4|^^)(^^|(i _ A)V4 



||(l-A)i/4^,||2 ||(l-A)l/4(^i||2 

(1 - A)V4^W (1 _ A)V4 (1 _ A)i/4|^^)^^^|(i _ A)i/4 



Tr(l-A)i/4^«(l-A)i/4 
Tr(l-A)V4^W(i_A)i/4_||(i_A)V4^^|| 



1 , 1/2 



which concludes the proof of the theorem. 



:i-A)l/4<^t||2 
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6 Control of the growth of the kinetic energy 

The goal of this section is to prove Proposition 15.2^ which gives control on the growth of the expec- 
tation of the kinetic energy operator /C with respect to the fluctuation dynamics UN{t; 0). 



Proof of Proposition \5.S\ Recall the definition (j4.5p of the unitary maps Ujsr{t; s) describing the evo- 
lution of the fluctuations; note that the generator C^it) of ZYiv(t; s) is defined in terms of the solution 
ift^ of the regularized Hartree equation 



idtfl 



(oiv) 



Vl^^ 



n 



("iv) 



1 



. + aN 



*l^l"-)pUl"-) 



(6.1) 



In the rest of this section, we will use the shorthand notation (pt = Vt ^ • By ()5.ip . by the assumption 
(/? G H'^{M.'^) on the initial data, and by Corollary 12.31 there exists v = v(T,k, \\^\\h'^) such that 



uniformly in N 



sup \\4>t\\H'^ < V 
\t\<T 

ifi'^ depends on N through the cutoff a^)- 



(6.2) 



We compare the growth of /C along the fluctuation dynamics Un and along a new dynamics Wat 
defined through the equation 



idtWN{t]s) = MN{t)WN{t]s) with >VAr(s;s) = l for all s G 



(6.3) 
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with the time-dependent generator 



M 



Nit) ■■= / dxa*.(l-A^)i/2a^. - X f dx {--^ * Ic^-ip) a*a^ 

-A / dxdy- — '4)t{x)(i)t{y)a*ax 

J \x-y\+aN " 

X f 1 — — 

-- dxdy- — {(j)t{x)<j)t{y) ala* + (ptix) cptiy) a.^ay} (6.4) 

z J \x — y\ + a]\f 

dxdy , _ , {My) o*xa*y I^nW) ax + Mv) «x '^^n{J^) a^ay] 



VN J \x-y\+aN 

X f 1 

-— : dxdy — ala* l^N{^f) aya. 

2N J \x - y\ + UN ^ 



y^x 



where l^Ar(s) is the characteristic function of the interval {—oo,'dN], that is l^]\f{s) = 1 if s < ^N 
and l^Ar(s) = otherwise. Here i? < 1 wih be fixed later to be sufficiently small. 

We split 

{UN{t; 0)n,}CliN{t; 0)f]> = (Wjv(i; 0)0, /C VVjv(i; 0)n) + {{UN{t; 0) - W^(t; 0))!^, /C VViv(t; 0)^) 

+ {UN{t; 0)n, )C{UN{t; 0) - W^(t; 0))J^) 

< (Wjv(i; 0)^, K, W^(t; 0)f]> + ||/C W7v(t; 0)(^|| || {UN{t] 0) - VV7v(t; 0))1]|| 

+ ||/CZYiv(i; 0)O|| II {UN{t] 0) - W^(t; 0))J7|| . 

(6.5) 

Proposition 15.21 now follows from Proposition 16. H Proposition 16. 2| and Proposition 16.31 and from the 
assumption that N^aj\f — >■ oo for some /3 > 0. D 

The first ingredient in the proof of Proposition 15. 21 is a bound for the growth of the kinetic energy 
/C and of its square w.r.t. the cutoffed evolution WAr(t; s). This is the content of the next Proposition, 
which will be proven in Section 16.11 



Proposition 6.1. Suppose that the assumptions of Proposition \5.S\ are satisfied (hut here the as- 
sumption N^aN — 7- cx) for some /3 > will not be used). Let the evolution WAr(t;s) he defined 
according to h6.3\) . with generator \6.4^ , and suppose that i} > is small enough. Then there exists 
C = C{'d,T,K,\\ip\\jj2) such that 

{WN{t; 0)n, JC^ nV(t; 0)f^) < C (6.6) 

for allt£R with \t\ < T. 

The second ingredient to prove Proposition 15.21 is a weak bound on the growth of the expectation 
of fC^ with respect to the dynamics UN{t;s); the proof of the following proposition is given in 
Section I 



Proposition 6.2. Suppose that the assumptions of Proposition \5.S\ are satisfied. Then there exists 
C = C{T,K,\\ip\\H2) such that 

/ /\r2\ 

{UNit;0)n,IC'^UN{t;0)n) <C( iV^ + ^-J (6.7) 

for allteR with \t\ < T. 
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Finally, we need to compare the two dynamics UN{t]s) and >VAf(t;s). The next proposition is 
shown in Section 16.31 



Proposition 6.3. Suppose that the assumptions of Proposition \5.2\ are satisfied. Let the evolution 
WAr(t;s) be defined according to 116. 3\) . with generator ^6.4^ , and with -i? > 0. Then, for any k £N, 
there exists C = C{k,-d,T, n, \\^\\h2) such that 



UMit;0)-WN{t;0))n\\<^ (1 + ^] (6.8) 



for allteR with \t\ < T. 



6.1 Growth of /C^ with respect to regularized dynamics 

The goal of this section is to prove Proposition 16.11 We will make systematic use of the bound (16.2 
(recall that in this section we use the shorthand notation (pt = fi )■ 



Observe that, by the definition (j6.4p of Mi^{t), we have 

+ ( / dxdy- — '4>t{x) (ptiy) a*aA 

\J \x-y\+aN " / 



+ {\ I'dxdy- \- {M^)My)<a; + h.c.}y (6.9) 

\2 J \x - y\ + ON ) 

+ \~rT^ / dxdy^ — {0t(y) a*a* l^Ar(7V) a^; +h.c.}) 

+ (wt:? dxdyr — ala*li)NiJ^)ayaA 

\2I\ J \x — y\ + On / 

where h.c. denotes the hermitian conjugate. First of all, we note that, by Lemma 16.41 below, the last 
term is bounded by 

(tttt dxdy- — ala*UN{^)aya^ < ^^/C^ 

\2N J \x — y\ + UN ' 

Therefore, choosing t? > sufficiently small, we find 

+ ( / dxdy I — 0((x)(/>i(2/)a*a^) 

^J \x-y\^aN " ) 

+ (- dxdy- — {(l)t{x)(j)t{y)ala*+h.c.}) 

\2 J \x — y\ + aN " -• / 

+ \~rTT / ^xdy- — [(t)t{y) ala* li)N{N) a^ +h.c.} 

V ^N J \x - y\ + aN ^ 

To bound the second term on the r.h.s. of the last equation, we note that 

d^ (^ *\'t>t?] 4ax<sup(^*|</.in M < \\(t)t\\]ii/2M . 

\\-\+OLN J X Vl-I / 
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(6.10) 



Since moreover A^ commutes with the operator on the l.h.s., we conclude that 



dx(^ *\<PtA ala,) <M\ (6.11) 



On 
Analogously, the third term on the r.h.s. of (IG.lOp is bounded by 

1 - \2 



f dxdy — 0i(x) My) «;a,) < Af' . (6.12) 

Next, the terms on the third line of (jG.lOp can be controlled as follows. Let 

A= dxdy — (j)tix)My)ala* . 

J \x-y\+aN ^ 

Since {A + A'f < 2{AA* + A* A), we find 

{'iIj,{A + A)^)< / dxdydx dy — — — {tp, a^a a^>ay'^) 

J \x-y\+aN \x' -y'\ + aN " 

, f A A A 'A ' Mx)My) '4>t{x')'4>t{y') , , ^ * * ^ ,\ 
j,^,\M^)Ut{y)\ \H^')Ut{y')\ 

a^v 1^;' — y'\ + oiN 
+ / dxd,dx' |^,(x)|^^M^L ^M^:L „«„ ^„ „a.. 



< ' dxdydx dy -. j -j— ; -. ||«xa«V'll llox'flv'^ 

\x — y\+ On \x — y \ + a^ 



X — y\ + ttN \x — x'\ + a^v " ^ 
+ / dxdy 



(|x-y| + a^)2ll^ll 

for arbitrary ^ G J^. Here we used that 

[a^Oy, al^a*,] = a*,ax5{y - x') + a%,ax5{y - y) + a*,ay6{x - x) + a*,aj,(5(x - y') 
+ 5{x - x')6{y - y') + 6{y' - x)5{y - x'). 

With Schwarz inequality, we obtain 



(V, {A + A*f^l,) < / dxdydx'dy' r 12 ll«-' V^ 



1*^2.,A < / A^A..A^'A..' \M^)nMy)\' 

{\x-y\+aNf "" ^ 



x-y|+aAr)2 "" J (Ix-yl + aN)' 

< 



^nMllfi||'At|r(V',(AA+l)» 



Thus 



.2 
Now, we estimate the terms on the fourth line of (|6.10p . Let 



If 1 

B = -j= / dxdy p- ala*(j)t{y) l^Af (AA)a2 

VN J \x-y\ + aN " 



3.13) 



i [dxdy- \- {</.t(x),/.i(y)<a; + h.c.})'< {M + if . (6.14) 

2 J \x - y\ + ON ^ ^ / 
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Then {B + B*f < BE* + B*B. The term BB* can be bounded by 

(V,i?i?*v> = ^ fdxdydx'dy'- \- — — — ]— — My)My') 

' N J \x - y\ + UN \x' -y'\+ aN 

X (l^Ar(A/' - 2) tp, ala*yaxal,a^'ay> l^Af (A/" - 2) tp) 

If 1 1 — 

= — / dxdydx'dy',- — — Mv) Mv') ,^,,, 

NJ \x-y\+aN\x' -y'\+aN (6.15) 

X (l^Ar(7V - 2) V, ala*yal,ax'axay/ l^Ar(7V - 2) ip) 

If 1 1 — 

+ — / dxdydy'i — — I — — My) My') 

NJ \x -y\+aN \x -y'\+aN 

X (l^7v(AA-2)V,a*a;a^ay 1^^(AA-2)V) 
for every ip G J^. From Schwarz inequahty, we find 

{^,BB*ij) < 1 fdxdydx'dy',--^^^ -^ ||a,.a,ay l^jv(AA - 2) Vf 

' NJ (\x-y\+aN)^ 

+ 4 f dxdydy' --^^^ ^||a,a,,l^^(AA-2)V^f 

NJ {\x-y\+aN)^ 

< 1 (sup/dy^iM!") [dxdx'dy'\\a,,a,ay,UN{Af-2)^^ (6.16) 

N \ X J \x-y\y J 

+ ^Lnv j dy^-^^^^ j dxdy'\\axay,UN{M -2)tljf 

The term B*B, on the other hand, is given by 

1 /" 1 1 — 

{ip,B*Bip) = — dxdydx'dy' — — — My)My') 

NJ \x -y\ + aN\x' -y'l+ON 

X (l^Ar(A/' - 1) ip, alaxayal,a*,ax' l^NiJ^ - 1) V') 

1 /■ 1 I — , 

= — / dxdydx'dy'-^ — — — Pt{y) Mv) ,^ _, 

NJ \x -y\ + aN\x' -y'\+aN (6.17) 

X (l^Ar(7V - 1) V', alal,a*y,axa^/ay ltfAr(7V - 1) V') 

1 /■ 1 1 — , 

+ — dxdydx'dy'.^ — — j— My) My') 

NJ \x -y\ +aN \x' -y'\+aN 

X (l^Ar(A/' - 1) V', < [axtty, al,a*y,] a^' ltfAr(AA - 1) ip) 
The first term on the r.h.s. is bounded in absolute value by 

— / dx dy dx'dy' .p^ r^-j- ry lla^ajj/a^' l^iv(A/' - 1) '0||||axax'aj, l,9Af(A/' - 1) V'll 

N J |x — y\ \x — y \ 

< ^ (sup I dy ^^) / dx dx' dy' ||a,a,my 1^^(AA - 1) V'f (6-18) 

<^Utfmm + l)'^'UN{M-l)ipf. 

When we insert (j6.13p in the second term on the r.h.s. of (j6.17p . we obtain contributions quartic in 
the creation and annihilation operators of the form 

If 1 1 - 

— / dxdydy'.^ — p^{y)My') {^^n{J^ - 'i-)'P,a*xa*,axayl^Ni-^f - 'i-)'P) 

NJ \x-y\+aN\y-y'\+aN y « 
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whose absolute value can be bounded by 



—7 dxdy dy' -. j -. -. ; 

NJ \x-y\+aN\y-¥\ + aN\ 



4't{y)My') i^mi-^ - 1)^, <ay OxOx' i^^a^ (aa - i)tp) 



If 11 

< — / dxdydy' —-—--— ——\(l)t{y)\ \(l)tiy')\ Way^a^l^Ni-^f - IjV'll hxax' 1^7v(7\A- 1)^\\ 

-'' J \^ y\ \y y \ 

\ y J \ y hf \ /J 

+ — f sup / dy— -^|0i(2/)n \ dxdy Hay/a^. I^tvCAA - l)?/^|p 

<^ll'At||^i||(AA + l)l^^(AA-l 



The other terms arising when we insert (j6.13p in the second summand on the r.h.s. of (j6.17p (both 
the quartic and the quadratic terms) can be bounded analogously. Together with (16.18p . we conclude 
that 



N 



From (j6.15p and (j6.19p . we find 



^Vn J 



dxdy ■ 



(y) <,< UN{^f) a. + h.c.} f <^{Af+lf UN{Af - 2) 



\x — y\ -\- ajv 
Combining (J6TT]) . (f6T2|) . ([6Tl]) . and (fOOj) we conclude (since I^at < 1) that 

IC^ < M%{t) + {M + if . 
Next, we observe that there exists a constant C = C{T, k, \\ip\\jj2) such that 

{WN{t; 0)n, {Af + if WN{t; 0)0) < C 



(6.19) 



3.20) 



3.21) 



(6.22) 



for all |t| < T. The proof of this bound is analogous to the proof of Proposition 14.31 (see Lemma 3.5, 
with M = tDN, and its proof in |21)). The only difference is that the generator M^it) contains a 
cutoff also in the quartic term (while in Proposition 14.31 the cutoff appeared only in the cubic term 
of the generator Ai]y{t)); this difference does not play any role in the proof of ()6.22p because the 
quartic term (with or without cutoff) commutes with the number of particles operator M (and thus 
with its powers). 

Finally, we control the growth of the expectation of A^^(t). To this end we compute, using (16. 3p . 



dt 
and thus 



Wiv(t; 0)n, M%{t)WN{t; 0)0 ) = ( W7v(t; 0)0, MN{t)MN{t) + MN{t) MN{t) >Vjv(i; 0)0 



^ (W^(t; 0)0, A^|,(t) Wiv(i; 0)0^ '^' 



< /Wiv(t;0)O,A^^(t)W7v(t;0)O 



3.23) 
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We have 



MNit) = - X dx 



1 



X dxdy ■ 



+ UN 
1 



* {<Pt4>t + <t)t4>t) 1 [x] ala^ 



|x — yj + ctAT 
A / dxdy 7- , , ((j)t{x)<t>t{y)a*yal + 'h.c. 



X 



X — y\ + a]\f 
1 



(6.24) 



dzdy 



(j)t{y)ala* l^Af (7\A) a^ + h.c. 



'N J " \x-y\+aN 

Observe that from the (regularized) Hartree equation ()6.ip we easily find that \\4>t\ 



< 



'^t\\H 



1 and 



|V</.i|| < \\{l-A)<Pt\\^ + 



V 



1 



+ a 



< 



< 



h\\m + 



* m 



k\\6 + 



1 



* (pt 



V^<Pt 



.25) 



mH^ + ll0tllf7i/2 \m\\H 



This implies, by (j6.2p . that there exists a constant C = C{T, k, W^Wh^) such that 

\\Mm<C for ah t G M with |t| < T. 



(6.26) 



Next, we bound the square of the terms on the r.h.s. of (|6.24p . Similarly to ()6.1ip and (|6.12p . 
we find 

1 



dx 



ttN 



* {4>t4>t + 4>t4>t) (x)a*.aa; ) < \m\\m\m\\H 



lA/-^ 



and 



dxdy ■ 



|x - y| + ON 
Moreover, similarly to ()6.14p . we obtain 



Mx)My) + Mx)My)] <Ajax ) < ||0t||//i||(/'t||/^iA/'^. 



(6.27) 



(6.28) 



dxdy ■ 



1 



</>t(x)(/>t(y) + <At(^)'/'t(y) aX + h.c. <Ut\\mUm+^f- (6-29) 



|x — y| + aN 

Finally, analogously to (I6.19P (replacing 4> with (j)) we have 
If., 1 



dxdy ■ 



{y)ala:, l^jv(AA) a, + h.c. < -\\Uh^ {N + 1) 



1 



'N J " \x - y\ + aN """" -^ ^ ,„,. .,-. ,--^ _^, 

From (fOiD . dnSZD, dnSHD, dnSSD, (IQUD . we find, using (f6:22D . that 

(W7v(t;0)f^,A^;vWW^(*;0)f^) < 1. 
Eq. (j6.23p then implies that there exists C = C{k,T, ||<^||//2) such that 

for alH G M with \t\ <T. Proposition O now follows from (jOT]) and (jO^j) . 



(6.30) 



D 
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Lemma 6.4. There exists a universal constant C > such that 

( [dxdy- ^-^ala;UN{M)aya,) < C ^^ IC" (6.31) 

for all a, -i? > 0. 

Proof. Denote 

V = / dxdy — — ala* l^Ar(A/')aya^ . 

J \x-y\+a " 

Then V (and thus V^) leaves the number of particles invariant and, on the n-particle sector, we have 
(^')^"^ = f^ E u._l.M J ' if-^^^ (6-32) 



N ^-^ \xi — X,- + a 

and (V^)(") = when n > -QN . Using the operator inequahty (see, for example, Lemma 9.1 in [2]) 

i-^ < (1 - A.)V2 (1 _ A,)V2 

we find 



\x — yp 



l<t<i<n V' ' ■?' ^ l<j<i<n 

2 7n" X2 " " (6-33) 

V 1=1 ^ 



and the lemma is proven. D 

6.2 Weak bounds on growth of K? with respect to fluctuation dynamics 

In this subsection, we show Proposition 16.21 Again, we will need the estimate (j6.2p : recall also that 
in this section we use the shorthand notation 0^ = (/jj ^ for the solution of (j6.ip . 

We write 
(Z^7v(i;0)J^,/C2^/^(t;0)rj) 

= / dxdy l^U^it- 0)f^, (1 - A,)V4 a* (1 - ^^fl^ a, 

X{1 - Ay)^/^ alii - /\yY'^ ayUN{t-^m) 

dxdy (n, (1 - A,.)VX(i; 0) 4 ^/jv(t; 0) (1 - A^fl^UUt; 0) a,. UN{t; 0) 
X (1 - Aj,)i/^^/^(t; 0) a*yUN{t- 0) (1 - A^)i/^Z^iV(*; 0) ay Unit; O)!^) . 
Next we use that (see (|4.7|) ) 

U*N{t;0)a^UN{t;0) = W*{^ip)e''^'^\a.^ - VNMx))e~''^''^W{VNip) 
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to conclude that 

-2 



X (1 - A,)V4 (a; - VnUv)) (1 - A,)^/^ {ay - Vn^v)) e-''^"N'W{VNip)n) . 

(6.34) 
For f £ L'^{R^), let 

vr(/) = a* (/) + a{f) = J dx (/(x) < + /(x) a,) . 

Then, from (|6.34p . we obtain 

{UN{t;0)n,IC^UN{t;0)n) 

+ 2\/iVRe (e-^^^* Pr(\/iV(/j)rj, /C 7r((l - A)!/^*) e-*^^*t^(\/iV(/p)0) 

+ 2A^ \\(t>t\\]ji/2 {e-'^N^ W{^/N^)VL, /C e-*^^*VF(\/]Vv3)S7> 

+ iV(e-^^^* VF(\/iVv9)f^, 7r2((l - A)^'^<Pt) e-''^'^«'W [Vn ^)n) 

Using Schwarz inequality, we find 

+ N{e-''^Ni W{^/Nip)^, 7r2((l - A)^/^*) e'^'^'^^^W {Vn ^)^) (6.35) 

+ iV'||0i|lW- 
Next, we observe that, for arbitrary / G L'^(R^) and ^ S J-", 

(V', vr2(/)V.) = ||7r(/)V'f < ||a(/)Vf + ||a*(/)Vf < ||/f (V", (AA + 1)V) • (6.36) 

Moreover, we have 

\ I 1 

%% = /C — V, where V = ^--r / dxdw ■; ; Qla*ayax ■ 

^ 2N J ^Ix-yl + a'^y^ 

Since 

~ Na 



and since [V,AA] = 0, we conclude that 



/c^ < {n^f + v'< {nPr + j;^M' (6.37) 



Inserting (|6.36p and ()6.37p in (|6.35p . we find 

{UN{t;o)n,ic'^UN{t]0)n) 



< {w{VNip)n, {n%f w{VN^)a) + -l-^(w{VN^)aM^ w{Vn^)q) 

+ NUt\\m/2 {W{^^)n, {M + 1) W{^^)n) + N^Ut\\Hi/2 (6.38) 

+ NUt\\m/2 {W{VN^)n, {M + 1) W{^^)n) + N^(t>t\\]j^„ . 
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< / \\{CN{s)-MNis))WN{s;0)n\\ds. 
Jo 



We recall that 



X f I 

MN{t)-CN{t) = —^=- dxdy- — (/)t(y) a*.a* (1 - l^Ar(A/')) a^ + li.c. 

X f 1 

+ ^7^ / dxdt/- ala*{l-l^Ni^f))aya^. 

2N J \x - y\ + ON '^ 



(6.40) 



Using the properties of Weyl operators listed in Lemma 13.21 it is simple to check that 

{W{VNip)Q,{M + l)W{VNip)Q) <N and {W{^ip)Q,M^W{VNip)n) < N'^ . (6.39) 
Moreover, we have 

= J dxdy {WiVNv)n, (1 - A,)i/4a* (1 - A,)^/^a, 

X (1 - Ay)^/^a*y (1 - Ay)^/^ay W{VNv)n) 
= j dxdy {n, (1 - A,.)'/'« - Vn^x)) (1 - A,)V4(a^ _ ^^{x)) 

X (1 - A,)i/4(a* - ^/iV(^(y)) (1 - Ay^/^ay - y/N^{y)) W{^ip)n) 

Inserting (j6.39p and (j6.40p into (j6.38p . and using the bound (j6.2p . we conclude that 

This completes the proof of Proposition 16.21 D 

6.3 Comparison of fluctuation dynamics with regularized dynamics 

In this section, we prove Proposition 16. 3[ 

We rewrite 

\\{UN{t]f))-WN{t;f)))n\\ = ||(l-Z^iV(i;0)W7v(t;0))J]|| 

ft _ _ 

< / \\UUs-M^N{s)-MN{s))WN{s-,^)^\\ds (g^4^^ 



(6.42) 



Analogously to (j6.20p . but with l^jy replaced by 1 — I^n, the square of the terms on the first line of 
last equation can be bounded by 

1 /■ 1 ^ x2 



/ dxdy-^ — (/>t(y)a*a* (1 - 1^N{-N'))ax + h.c, 

J \x-y\+aN " 



N J \-^ — y\-r uLi\i / /g ^„N 

<l(AA + l)3(l-l^^(AA-2)). 
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As for the second term on the r.h.s. of (j6.42p . its square can be estimated as follows. 



If 1 

— dxdy — ala* (1 - l^Ar(7V)) Uya^ 

NJ \x-y\+aN ^ 

= 1^ / dx dy dx'dy' — — — 

N^ J \x -y\ + aN \x' - y'\ + un 

X (1 - ltfAf(AA - 2))ala*ayaxal,a*y,ay:a^>{l - l^Ar(A/' - 2)) . 
From ayaxal,a*, = 0*^,0*, ayax + [0^0^,0*/ a*,], and from (|6.13p . we conclude that 



3.44) 



— - dxdy j 

NJ \x-y\+aN 



ala* (1 - l^Ar(7V)) Qytta 



1 



< 

- ATS 



dx dy dx'dy' 



1 



+ 



+ 



1 



X - y\ + UN \x' - y'\ + UN 
X (1 - l^Ni^f - 2))alayal,a*y,ayiax'ayaxil - l^dNi-^f - 2)) 

dx dy dy' 



N"^ j " " \x — y\^ a^\x — y'\^ afq 

X (1 - \^n{M - 2)) alaya*>ayiayax (1 - 1^n{J^ - 2)) 



(6.45) 



1 



Ar2 



dxdy 



1 



(|x-y| + aNy 



(1 - UN{^f - 2)) alalayttx (1 - 1^7v(AA - 2)) 



<^pr^{M+l)'{l-UN{J^-2)) 



iV2 



a 



N 



From ([03]) and (j05]) . we find that 



{MN{t)-CN{t)f<ii+ ^ 



N N^a% 



(AA+l)^(l-l^jv(AA-2)) 



.46) 



For every /c G N, we have (1 - Un{M - 2)) < {M - 2f/{'dNf. Therefore 



{MN{t) - CN{t)y 



< 



+ 



1 \{M+l 



^ ' N ' ma^ 



N 



\k+i 



{§Nf 



(6.47) 



Analogously to Proposition 14.31 (see Lemma 3.8 and its proof in [21j), there is C = C{k, k, T, \\(p\\h'^] 
such that 

(Wjv(i; 0)^, (AA + 1)^+^= Wjv(t; 0)^) < C 



for all \t\ < T. From (j6.46p . we find that, for every A; G N, there exists C = C{-d, k,T, k, \\^\\h^) such 
that 

{MN{t) - /:N{t))WN{t;0)n - ^ '^ 1 



- N^\N^ Nun) ' 



The proposition now follows from (j6.4ip . 
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